Observation of Antinormally Ordered Hanbury-Brown— Twiss Correlations 
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We have measured antinormally ordered Hanbury-Brown— Twiss correlations for coherent states 
of electromagnetic field by using stimulated parametric down-conversion process. Photons were 
detected by stimulated emission, rather than by absorption, so that the detection responded not 
only to actual photons but also to zero-point fluctuations via spontaneous emission. The observed 
correlations were distinct from normally ordered ones as they showed excess positive correlations, 
i.e., photon bunching effects, which arose from the thermal nature of zero-point fluctuations. 
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Since Planck's quantization hypothesis of electromag- 
netic field and Einstein's photoelectric theory appeared, 
the quantum nature of electromagnetic field has been in- 
tensively studied 0. The normally ordered photodetec- 
tion theory of Glauber Q played a central role in these 
studies [lj . This theory provided a formal explanation of 
the Hanbury-Brown— Twiss (HBT) type correlation mea- 
surements |3j, in which the electromagnetic fields were 
detected at two separated space-time points. The nor- 
mal ordering reflects the fact that the electromagnetic 
field is detected by an absorption process. Thus, such 
a photodetection is insensitive to zero-point fluctuations 
because the photodetection probability for the vacuum 
state is zero, i.e., (0|a^a|0) = 0. This insensitivity is 
the very reason why Planck's spectrum of black-body 
radiation is always convergent, even if the electromag- 
netic energy caused by zero-point fluctuations is diver- 
gent y|. It is also the reason why pieces of information 
on the field, more precisely, the vacuum components of 
the field's density matrix, are lost during the detection 
process and thus the initial density matrix cannot be log- 
ically reversible 5], that is, cannot be calculated from the 
post-detection density matrix and the readout of the de- 
tection. 

When photons are detected by stimulated emission, 
however, antinormally ordered photodetection can be re- 
alized [a|. Since the detection responds not only to 
actual photons but also to zero-point fluctuations via 
spontaneous emissions, the system's information can in 
principle be conserved during this photodetection pro- 
cess • The photon-counting statistics of the detection is 
then distinct from that of the standard normally ordered 
photodetection especially in the region where the aver- 
age photon-number in the concerned modes is small 0. 
Although the emission-based antinormally ordered pho- 
todetector was originally proposed by Blocmbergen, as 
the quantum counter, as early as 1959 Q, no report on 
the experimental realization exists till date to the best of 



our knowledge. 

As a demonstration of the antinormally ordered pho- 
todetection, we realized the quantum counters by using 
stimulated parametric down-conversion process and mea- 
sured the HBT correlations for coherent states in antinor- 
mal order. Here, we utilized an ultrashort pulsed laser as 
a pump field for the parametric process to obtain large 
nonlinear response of the crystal and to overcome the 
slow response time of the detectors 0, . The observed 
correlations deviated from the standard normally ordered 
ones as they exhibited excess positive correlations, i.e., 
photon bunching effects. The deviation can be attributed 
to zero-point fluctuations, which are known to possess 
thermal characteristics [tj [lfj . 

To clarify how to measure the antinormally ordered 
HBT correlation, we pursue the time evolution of a 
single- mode annihilation operator, a\ n , as the signal 
field. A schematic illustration of the correlator is 
shown in Fig. 2] (a). First, operator a\ n is coupled 
with an operator, bi„, via parametric interaction with 
the pump field and evolved into a "ut = <Hn cosh[sL] — 

bi n e" 5 sinh[sL]; while operator bi n becomes b out = 
bi n coshfsi] — aln e 1 "® sinh[s_L] under a perfectly phase- 
matched condition Here, L is the length of the crys- 
tal, and several parameters in the interaction, such as the 
intensity of the pump field and the second-order non- 
linear susceptibility of the crystal are included in pa- 
rameters s and d. An energy diagram of this process 
is shown in the inset of Fig. (a). Next, by divid- 
ing the field represented by b ou t into fields b\ and 62 
with a half-wave plate and a polarizing beam-splitter, 
the standard HBT interferometer for field b out is formed 
as shown in Fig. (a). Here, the imperfect quan- 
tum efficiencies of the photodetectors and several op- 
tical losses are taken into account and modeled by in- 
troducing auxiliary beam-splitters with vacuum fields v\ 
and v 2 l\- Output fields d\ and di of Fig. (a) are 
then written as d\ — ^/rji (T b ou t +H v) + i\fl — r/i vi 
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FIG. 1: (a) Schematic illustration of the antinormally ordered HBT correlator based on stimulated parametric down-conversion, 
(b) Experimental setup (the mode-locked Ti:Sapphire laser and the frequency doubling crystal are not shown but placed at the 
left side of the figure). The correlations for the signal field, din, are acquired by counting the delayed- coincidental photodetection 
events of fields b\ and 62 with detectors 1 and 2. 



and di = yjryi (TZ b out + T v) + - 772 U2, respec- 
tively. Here, 771 and 772 are the total photodetection 
efficiencies for detectors 1 and 2, respectively; TZ and 
T denote the reflection and transmission coefficients at 
the beam-splitter, respectively, and can be varied with 
the half-wave plate and the polarizing beam-splitter; 
and v is an auxiliary vacuum field introduced from the 
empty port of the polarizing beam-splitter. Since the 
modes relevant to operators bi n , V, vi, and V2 are ini- 
tially vacua, the surviving contribution to the number 
of the normally ordered photodetection events for fields 
di and d 2 coincides with that of the antinormally or- 
dered photodetection for field aln up to a constant fac- 



tor, i.e., (ndj 
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and (n d2 ) = (d 2 d 2 ) = 77 2 |72.| 2 sinh [sL](d} n dl n '), re- 
spectively. Here, the angle brackets indicate quantum- 
mechanical expectation values. Furthermore, the num- 
ber of coincidental photodetection events of fields 

di and ri 2 results in {hd 1 hd 2 ) = (d\ d\d 2 d 2 ) — 
ViV2\1~\ 2 \7l\ 2 sinh 4 [s L](al„a' in al n ^ ' a ir J) , where we use the 
commutation relation for each operator and relations 
T1Z* = -TIT* and cosh 2 [.s] - 1 = sinh 2 [s]. Thus, the 
surviving contribution to the coincidences turns out to 
be made only by operators aln and dl n in antinormal 
order. Consequently, we can evaluate an antinormally 
ordered HBT correlation for field a\ n as follows: 
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Note that Eq. holds regardless of the splitting ratio 
at the beam-splitter, the quantum efficiencies of the de- 
tectors, and the optical losses. 

To analyze the correlation more realistically, the de- 
tection process should be treated with time-dependent 
and continuous- mode field operators Q. In this treat- 
ment, the antinormally ordered HBT correlation, Eq. 0J, 
becomes time dependent. Moreover, since the response- 
time jitter of the detector is larger than the pump pulse 



duration (the duration of the parametric interaction) but 
smaller than the time interval between two successive 
pulses, the relevant information on the time dependence 
is embodied in the integrated number of delayed coinci- 
dences over the response-time jitter of the detectors 0, El ■ 
Thus, what we should measure to evaluate the antinor- 
mally ordered HBT correlation is the correlation of the 
pulses of mth neighbor: 
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where t and t + mr are the initial time of the integra- 
tion for detectors 1 and 2, respectively, T is the duration 
of integration, and r corresponds to the time interval 
between two successive pump pulses. Here, since the co- 
herence time of the down-converted field is far shorter 

= 1 for 
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Note 



than the pump pulse interval, we have g 
m 7^ 0, thus we only focus on the value, g 
that the duration, T, can be approximated by infinity be- 
cause the response-time jitter of the detector is far longer 
than any other relevant time scale. By taking account of 
the pump spectrum and the bandwidth of the interfer- 
ence filter [tJEl, the correlation, g^ A ^ of Eq. for a 
coherent state becomes 12] 
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(h) + i [ (n) + 1 ; 



(3) 



where (h) = (<zj n '<4 n ) is the average photon number of 
the signal field. The parameter, a, represents the indis- 
tinguishability of two emitted photons, which are respon- 
sible for a coincidental count. Here, the value, a = 1, 
corresponds to the case that two emitted photons are 
completely indistinguishable. 

A rough sketch of our experimental setup is shown in 
Fig-ffl(b). A parametric down-converter was formed with 
a beta-barium borate (BBO) crystal (thickness: 2 mm) 
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FIG. 2: Experimental results, (a)-(c): The number of sin- 
gle counts in detector 1 as a function of the optical path- 
length difference between the pump field and the input signal 
field, ain, whose average photon number (n) is (a) (vac- 
uum), (b) L09, and (c) 7.72. (d)-(f): The number of delayed- 
coincidental counts in detectors 1 and 2. Here, the coherent 
field with the average photon number indicated in each of the 
figure on the left was used as input signal. The insets indicate 
the accumulated number of coincidences within 3 ns. 



and a pulsed pump field (wavelength: 399 ran, average 
power: 195 mW, pulse duration: 100 fs, and repetition 
rate: 82 MHz) from the second harmonic of a mode-locked 
TkSapphire laser (Spectra-Physics Tsunami). The crys- 
tal was arranged to be type-I phase-matched so that ver- 
tically polarized frequency-degenerate down-converted 
photons (wavelength: 798 nm) were spontaneously emit- 
ted as the spatially- nondegenerate fields, ciout and b ou t, 
separated by ±2.54° with respect to the pump. A heav- 
ily attenuated coherent state (a laser output with the 
fundamental wavelength (798 nm)) was used as a sig- 
nal field represented by ain. Its average photon num- 
ber can be varied by the attenuator in Fig. ^ (b) . The 
signal field, <xi n , was amplified via stimulated paramet- 
ric down-conversion when the optical paths of the sig- 
nal and the pump fields were properly adjusted with 
the retroreflector and the mirrors shown in Fig. ^ (b). 
One of the stimulated fields, b out , was further split into 
fields bi and 62 with a half-wave plate and a polarizing 
beam-splitter as shown in Fig. 2] (b). After propagat- 
ing about 800 mm from the crystal and passing through 
an interference filter (bandwidth: 5.0 nm) each field was 



coupled into a single-mode fiber, which acted as a spa- 
tial filter, and then detected by Si-avalanche photodi- 
odes (Perkin Elmer SPCM-AQR-14; detection efficiency: 
55 %, response-time jitter: 350 ps, and dark-count rate: 
100 s _1 ). Electric pulses produced by the detectors were 
processed by digital- logic gates (Kaizu Works KN470), 
and the number of coincidences (threefold coincidences) 
in two (three) detectors as well as single counts in each 
detector were counted with pulse counters (Stanford Re- 
search System SR620) for aligning the setup. The pulses 
from detectors 1 and 2 were also fed into a time-interval 
analyzer (Yokogawa TA520; time resolution: 25 ps) for 
measuring their delayed coincidences. 

Figures |2Ia)-(c) show the number of single counts in 
detector 1 (during 1 s) as a function of the optical path- 
length difference between the pump field and the signal 
field, din, where average photon numbers (n) of the sig- 
nal fields for (a), (b), and (c), differ each other. Here, we 
set the half-wave plate in Fig. ^b) to be \T\ 2 — 1 and 
\K\ 2 = 0. The counts were enhanced by the stimulated 
emission within the area where the optical path-length 
difference was small, i.e., two pulses overlapped. On the 
other hand, the constant background counts, which were 
independent of the optical path-length difference, were 
attributed to the spontaneous emissions. The maximally 
enhanced counts due to the stimulations are given by 
{hdi) — rj\\T\ 2 sinh 2 [s]((n) + 1), whereas the counts with- 
out any stimulations are given by {hdx) — ili\T\ 2 sinh 2 [s]. 
Thus, by comparing these counts the signal field's aver- 
age photon number, (fi), can be quantitatively evaluated 
as (a) (vacuum), (b) 1.09, and (c) 7.72, respectively. 

Figures Eld)-(f) show the number of coincidental 
counts in detectors 1 and 2 (during 1800 s) recorded for 
various time delays (25-ps time bin). Here, the coher- 
ent field with the average photon number indicated in 
each of the figure on the left was used as input signal, 
and the half-wave plate was set to be |T| 2 = 1/2 and 
\1Z\ 2 — 1/2. The peaks of the coincidences at 12.2-ns in- 
tervals corresponded to the 82-MHz repetition rate of the 
mode-locked laser, and the width of each peak was dom- 
inantly determined by the 350-ps response-time jitter of 
the detectors. As mentioned before, the relevant infor- 
mation on the correlation should be extracted after ac- 
cumulating the number of coincidences within the width 
of each peak. The insets in Figs. I3d)-(f) indicate the 
accumulated number of coincidences within 3 ns, which 
is sufficiently larger than the response-time jitter of the 
detector. The second peaks in Figs. |2Jd)-(f) correspond 
to the coincidences at zero time delay (m = 0), where 
two photons produced by the same pump pulse are re- 
sponsible for the coincidences, and they are larger than 
any other peak. 

The antinormally ordered HBT correlations, g^ A ^ , 
were calculated as follows. As mentioned before, un- 
less the delay mr in Eq. J2J is zero, g^ A ^ can be re- 
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FIG. 3: Normalized antinormally ordered HBT-type corre- 
lation g^ A ^ ■ The open circles indicate the experimentally 
evaluated correlations with error bars for both axes. The 
dotted-line curve is the theoretical prediction obtained by the 
single-mode calculation. The full-line curve was obtained by a 
X 2 fitting of the data to the prediction with continuous-mode 
treatment. 



garded as one, i.e, no correlation survives. Therefore, 
the accumulated number of coincidences in the second 
peak (m = 0) normalized by that of the non-zero de- 
lay peak (m ^ 0) can be considered as the correlation, 
g&\. A Y> _ t evaluate the variance as well, we used five 
peaks (m = —1, 1, 2, 3, 4) at non-zero delay to determine 
the normalization factor. The averages of the correla- 
tions for seven input fields (the average photon numbers 
are 0, 0.71, 1.09, 2.29, 4.36, 7.72, and 10.61) are shown in 
Fig. as the open circles with vertical error bars (three 
standard deviations). Each horizontal error bar (three 
standard deviations) was evaluated by assuming that the 
number of photodetections followed the Poissonian dis- 
tribution. As opposed to the normally ordered HBT cor- 
relations for coherent states, the results show the photon 
bunching effects (g^ A ^ > ff^' = 1), which were trig- 
gered by the influence of zero-point fluctuations [9! Hof. 
The bunching effect becomes inconspicuous as the aver- 
age photon number of the input field increases and ap- 
proaches the normally ordered value, g — 1 |6|. In this 
sense, this discrepancy can be viewed as a purely quan- 
tum effect. From another point of view, these bunch- 
ing effects can be ascribed to the fact that a coherent 
state is not an eigenstate of a creation operator, i.e., a 
measurement operator of the antinormally ordered pho- 
todetection. Thus, the coherent properties of coherent 
states no longer hold in the antinormally ordered HBT 
correlations. 

The dotted-line curve in Fig. |21 indicates value g^ A ^ 
of Eq. (|3J) with a = 1, that is, the simple prediction 
obtained by Eq. |Q, where all the relevant fields are 



treated as single-mode states. The full-line curve in 
Fig. indicates the value of g 2 } with a — 0.45, which 
is obtained by applying a \ 2 fitting with a as a fit- 
ting parameter, i.e., by finding the minimum value of 

■ = Eli s|> - [i + "(t^ttt + itJM) 2 - Here ' 

(h)i, Hi, and a Vi are the average photon number, the cor- 
relation, and the standard deviation of the correlation 
for ith experimental value in Fig. EI respectively, and the 
variances of the horizontal axis are neglected. To take a 
quantitative look at how good the fitting is, we executed 
a \ 2 test of goodness-of-fit. The statistical distribution 
of value s is supposed to obey the x 2 distribution with six 
degrees of freedom (seven data — one unknown parame- 
ter), which has value x 2 (6) = 12.6 at the upper 5% point. 
On the other hand, with the best fitting (a = 0.45), the 
value, s, reaches 10.2, which is well below 12.6, thus, we 
can conclude that the experimental results are in good 
agreement with the continuous-mode analysis of the de- 
tection process. A detailed analysis shows that the 
resultant reduction of correlation, i.e., a — 0.45, can be 
mainly attributed to the imperfect erasing of the time- 
stamp information of the down-conversion processes (re- 
duces a to about 0.9), and the spatial-mode mismatching 
(reduces a further to about 0.55). 

In conclusion, to demonstrate antinormally ordered 
photodetection, the HBT correlations for coherent states 
have been measured in antinormal order by using stim- 
ulated parametric down-conversion. Since the measure- 
ment operator was no longer an annihilation operator, 
but rather a creation operator, even the coherent states 
exhibited the photon bunching effects. The emission- 
based antinormally ordered photodetection may provide 
an interesting alternative for monitoring quantum sys- 
tems owing to its sensitivity to zero-point fluctuations. 
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